Tunneling between Two Helical Superconductors via Majorana Edge Channels 
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We discuss electric transport through a point contact which bridges Majorana fermion modes 
appearing at edges of two helical superconductors. The contents focus on effects of interference and 
interaction unique to the Majorana fermions and role of spin-orbit interaction (SOI). Besides the 
Josephson current, the quasi-particle conductance depends sensitively on phase difference and rela- 
tive helicity between the two superconductors. The interaction among the Majorana fermions causes 
the power-law temperature dependences of conductance for various tunneling channels. Especially, 
in the presence of SOI, the conductance always increases as the temperature is lowered. 
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Fractionalization of electrons attracts recent intensive 
interest. A chiral fermion at an edge of Quantum Hall 
system is an example of fractionalized electrons. A fully 
gapped bulk state spatially separates right-going and left- 
going chiral fermion, which leads to the absence the back- 
ward scattering [l| . The robustness against disturbances 
such as disorder and interaction is a common feature 
of the fractionalized states. Therefore the fractionalized 
states are expected to have dissipationless feature which 
is a key property on application to quantum information 
processes Majorana fermion (MF) is another exam- 
ple of the fractionalized electron and has been recently 
discussed in the context of condensed matter physics [31 . 
Since its field operator in real space satisfies a relation 
7 = 7^, MF is often called real fermion and is regarded 
as a half (fraction) of a usual complex fermion. 

Superconductors and superfiuids are the most promis- 
ing candidates which host MF because the particle num- 
ber is not a good quantum number in these systems as 
required by the MF field. Actually the existence of MF 
has been discussed in a vortex core or at an edge of the 
chiral p -\- ip superconductor [3], ^He and Bose-Einstein 
condensates 0, Q, at an interface between a supercon- 
ductor and a topological insulator ML and in a quan- 
tum Hall edge (state) with = 5/2 Q. Chiral fermion 
modes appear only when the time-reversal symmetry T 
is broken [l|. Under preserving T-symmetry, the part- 
ner with the opposite chirality always coexists. In this 
case, the edge channels are referred to as helical. Along 
the edge of the two-dimensional quantum spin Hall sys- 
tems, the helical fermions appear [9|. Analogously, at 
the edge of the helical superconductors, the helical Ma- 
jorana fermions appear as the Andreev bound states fl^- 
|l2|. Noncentrosymmetric superconductors with domi- 
nant spin-triplet p-wave order parameter are a realistic 
playground of helical Majorana fermions [l^, [HI . In ad- 
dition, MF excitations are expected in topological super- 
conductors [13] • Thus understanding of novel phenomena 
peculiar to MF is highly desired [ij] . Although all of re- 



cent developments have assumed non-interacting MF's, 
effects of interaction among MF's have been an impor- 
tant open question. 

In this paper, we study theoretically low energy elec- 
tric transport through the MF modes appearing at edges 
of two helical superconductors with taking into account 
the interaction among MF's. The model of four in- 
teracting MF modes can be mapped into the spinless 
Tomonaga-Luttinger model by introducing two fictitious 
chiral fermions. This enables us to analyze low energy 
physical phenomena of a MF by using the bosonization 
technique. We show that the conductance depends sen- 
sitively on the relative helicity of two helical supercon- 
ductors, the phase difference of two superconductors, and 
the spin-orbit coupling at the point contact. These are 
the features of the interference effect unique to Majorana 
Josephson junctions. Our main results are summarized 
in Eqs. (HH), ^ and ([26|. 
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FIG. 1: (Color online) Schematic picture of helical Majorana 
excitations at edges of two helical superconductors. 

We consider the interacting helical edge channels fl^ 
as shown in Fig.l. By using solutions of the Bogoliubov- 
de Gennes equation, Majorana fermions at the edge of 
the helical superconductors are described by [12] 

Hq = - iv I dx [-fRjix)dx'yiij{x) 

~lLjix)d^-fLj{x)] (1) 
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where 7^(a;) for /z = {Rl, LI, R2, L2) are the four 
species of Majorana fermion field satisfying the anti- 
commutation relation {7^1(2;), 7/j' (a;')} ~ (l/2)(5^^^/5(a; — 
x'). The electron operator is expressed in the low energy 
sector as 

=e*'^^/'e*'^'x7,,a(a;), (2) 

r e"/4 (j,a) = (l,t)and(2,;) 
^ /4 (j, a) -(!,;) and (2, t), 

where (^j is the phase of superconducting order parame- 
ter for the two superconductors indicated by j = 1 and 
2, a =tii represents pseudospin index, and ipa- is the 
relative phase of the pair potential for the two pseu- 
dospins. In Eq. we have considered pair potential 
in two-dimensional ^He-B phase as an example of the 
helical edge state. As shown in Fig. [U the pseudospin 
index a —'[, J, is related to the moving direction of the 
chiral Majorana fermions {R,L), which are basically de- 
termined by the pair potential and the spin-orbit inter- 
action (SOI) in bulk region. 

In addition to Hq, we consider the two terms in the 
Hamiltonian. At first, the interaction Hint, comes from 
the electron-electron interaction as given by -ffci-int. — 
JdxJdx'Cl{x)Cl{x')V{x - x')Cp{x')Cc,{x)/2, where 
^a{i3) is the electron operator with a and /? labeling the 
electron spin. The original electron spin a is expressed 
by the linear combination of pseudospin a in the presence 
of the SOI. From the fact that ^'f^.^'io- cx (7io-)^ = const., 
and assuming that the overlap of the wavefunctions be- 
tween the two edges is negligible, the only remaining in- 
teractions derived from i?oi-int. in low energy are 

Hi = Ui J dx[^\^^ 1,1^1^^ 2A + h.c.] 

H2 = U2J da;[*l^^*i,^*2,t*2.i + /i-c-]. (4) 

Note here that more than four species of Majorana 
fermions and SOI are indispensable to having effective 
interaction [T6| and that interaction terms including spa- 
tial derivatives are irrelevant in the low energy limit. Ex- 
pressing Eqs. (U) by Eq. ([3]), we obtain 



^int. — . 



dx-fBlix)-fR2{x)jL2ix)-fLlix), (5) 



with g = -2Ui cos{2(ps) + 2U2- 

The last term is the tunneling Hamiltonian between 
the two edges represented by 



i/T - - to ^ [*L(0) {ao -I- l\ ■ <t}^^^, ^2,a'{Q) 



(6) 



+ *2,.(0){^o-*A-^},,,,vi,i^,,(0) 

=2ita [cos{ip/2)A- + A3 sui{ip/2)A+ 
- cos((^/2)A_B+ + sin((p/2)A+B_] , 
A± =7i,t72,t ± 7i472,4.> B± = 7i,t724 ± 7i472,t) (7) 



with A_ = Ai cos — A2 sin Lps and A+ = Ai sin Lp^ + 
A2COS(y9s, where tp = ip^ — ip2 is the macroscopic phase 
difference, Lps = 'f^ — f\. is the difference in the spin- 
dependent phases, cr = (cri, (T2, fs) are the Pauli matri- 
ces, and CTo is the 2x2 unit matrix. The width and the 
length of a point contact is indicated by a. In Eq. ([6|). 
we consider the SOI at the point contact described by 
A = (Ai, A2, A3). By applying electric fields at the point 
contact, it is possible to induce the Rashba-type SOI as 
A = gso{—Ez,Q,Ex), where gso is a coupling constant. 
Ex and Ez correspond to the potential gradient in the 
X and z direction, respectively. When the Dresselhous- 
type SOI can be introduced, A2 also becomes non-zero 
value. The operator of the electric current is calculated 
from the equation J = edt / dx^\ ^{x)'^ i^„{x). We 
find that a relation J{{p) = eHT{<f + T^) always holds. We 
assume that the pseudospin of the right (left)-mover is f 
(l) at the edge 1 as shown in Fig.[TJ Therefore we choose 

71. t — 7fli s-nd 71^4, = 7li. At the edge 2, we choose 

72. t = 7L2 and 72,^ — 7_r,2 for the equal helicity config- 
uration, and 72, f = 7fl2 and 72,^ ~ 7l2 for the opposite 
helicity configuration. 

To analyze the Hamiltonian, we introduce a com- 
plex fermion field by 'ipAix) — jai{x) + 17^2(2;) and 
iIj\{x) = 7^11(2;) — 17^2(2;) with A = R and L. These 
operators satisfy the usual fermion anti-commutation re- 
lations: {^ljA{x),i'A'{x')} = and {^pA{x),ipl^,{x')} — 
Sa,A'S{x — x'). We rewrite the Hamiltonian Hq + Hint, in 
terms of these complex fermion operators as 



Ho+Hint. = -iv 
9,A 



dx 



+ -^'^Rix)'^Ri^)^Li^)^Lix) 



V'fl(a;)(9xV'-R(a;) - V'l(a;)9a;V'L(2;) 

const.. (8) 



This Hamiltonian is exactly that of the massless 
Tomonaga-Luttinger model. It is extremely simple as 
compared to that of interacting helical edge fermion sys- 
tems [iTj . The combining two Majorana fermions corre- 
sponds to a spinless chiral fermion. Thus a chiral Majo- 
rana fermion is considered as a quarter fraction of a usual 
spinless fermion. By applying the standard bosonization 
technique [i8l[ll|, the complex fermion fields are repre- 
sented by boson fields as. 



d 

^<Pl{R){x) =: V'I(fl,)(a;)V'L(i?)(a^) (9) 
(l){x) = 0fl(x) + (t)L{x), 9{x) = (t)R{x) - (I)l{x), (10) 



V'i?(2^) = 



\J2-KaQ 



exp 
exp 



-{0(x)-K^(a;)} 
l{_0(^)+^(^)} 



(11) 
(12) 



where r\]i and 77^ are the Klein factor. Eq. ([5]) is then 



3 



transformed into 

Ho + H,,,,. = ^ fdx^^^^P^+K{d,0{x)}\ (13) 



V 1 



K 



' l-g/iSnv) 
l+5/(8H' 



(14) 



We first discuss the tunneling efi^ect in the equal helicity 
configuration. By using the bosonization technique, the 
tunnehng Hamihonian becomes 



+ 



ta 

TT 



A, 



■ sm 



dx(t>{x) 



x=0 



Trao 







cos 






.2. 



sin 0(0) - A3 sin | sin(/i(0) . (15) 



In Eq. (IT5|) . the terms including dxO{x) and dx(t>{x) 
represent the forward tunneling process: hopping of 
left(right)-mover to left(right)-mover. On the other 
hand, the terms including sin6'(0) and sin^(O) represent 
the backward tunneling: hopping of left(right)-mover to 
right(left)-mover. Before turning into the conductance, 
the Josephson current should be clarified. Within the 
second order perturbation expansion, we find 



J =eA 



at 


2 




sin (p 


TTV 



1 



A3 - -^K 



(16) 



where we have used (?y_R?7L)^ = —1, (ao)^^ = kp = A/u, 
and A is the amplitude of pair potential at sufficiently 
low temperature T ^ Tc with Tc being superconducting 
transition temperature. The ground state of junction is 
at = for A = [13 ■ Eq. does not recover a 
usual relation J cx (1 + A ) expected in s-wave Josephson 
junction even in the absence of interaction, (i.e., K — 
1). This is a characteristic feature of Majorana fermion 
excitation. 

On the basis of the linear response theory, we cal- 
culate DC conductance of the point contact using 
the standard Kubo formula, a = —limi^^Q+[Q^{Lo) — 
Q^{0)]/ {iuj), where the correlation function is obtained 
by Q^{^^) = Q{i^^n u! + iS) with Q(wn) — 

— /q'^^"^ dre*"""^ {J{t)J{0)), t is the imaginary time, and 
Lu„ is the Matsubara frequency. The following four terms 
contribute to Q{LOn), 

(J(r) J(0)) = {F^dx9iT)dxe{0) + F^dxHr)dxm),=o 
+ (Bj sme{T) sin 0(0) + sin ^(t) sin m) ' (1^) 

where Fg = cos((^/2)A+/2, F^ = sin((y3/2)A_/2, Bg = 
sin((p/2)/ao and B^ = cos{ip/ 2) X3 / ao . After deriving 
the effective action for 9{0, r) and (0, r), w e obtain the 
scaling equations by following Refs. |2l| - [23| . 



dl 
dBlil) 
dl 



.2 1 



1 



=2{l-K)Bl{l), 



(18) 
(19) 







A = 


A / 


Equal helicity 










K = 1 





const. 




K <1 





rp2K — 2 




K > I 





const. 


</p/0 


K ^ I 


const. 


const. 




K < 1 




rp2K — 2 




K > 1 


rj-i2/ K — 2 


rp2 / K — 2 


Opposite helicity 








(p = 


K = 1 





const. 




K < 1 





const. 




K > 1 





rp2/K-2 




K = 1 


const. 


const. 




K < 1 


const. 


rp2K-2 




K > 1 


const. 


rp2/K-2 



TABLE I: The temperature dependence of the most domi- 
nant terms in conductances given in Eqs. (|23p and (|26|) at 
low temperature, where ip is the phase difference between the 
two helical superconductors, A represents spin orbit interac- 
tion at a point contact, and 'const.' means the conductance 
independent of temperature. Here K = l{g = 0) represents 
no interacting case, while K < l{g > 0) and K > l{g < 0) 
are the interacting cases. The conductance depends on the 
relative helicity of the two superconductors: the equal helic- 
ity configuration (upper column) and the opposite one (lower 
column) . 



where fi = Ae~' gives a boundary between high- and 
low-frequency and A is the high-frequency cut-off. It is 
concluded that Be and are relevant for X > 1 (g < 0) 
and K < 1 (g > 0), respectively. 

The forward tunneling terms are calculated to be 



X—Xq / 



1 



X{r) 



K 

{dx4>{x,T)dx4>{xMx=xo) = kx{t), 
^W = ^^(.)-^T^, 



(20) 
(21) 

(22) 



Finally we reach the DC conductance for the equal helic- 
ity configuration 



a 

Gq 



^+ 2 

TT— cos 

K 



(I) .»(!).. 

nXlKsin^ (|) -|- A^ cos^ 



T 
T 



2/K-2 



2K-2 



(23) 



where Go 

rl/To 



= itaej'KV)^ and 

Ag/(,'^'^" drT(sinyl(r)sinyl(0)) ioT A ^ B and 
are the correlation function at T = Tq < Tc. Remark- 
ably, the conductance depends on the phase difference 
between the two superconductors. This stems from a 
peculiar feature of Majorana fermion. The two terms in 
the tunneling Hamiltonian, which are Hermite conjugate 
to each other, turn into the same form due to the "real" 
nature of Majorana fermion giving the interference 
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effect. At (p = 0, the conductance vanishes in the 
absence of the SOI, (i.e., A = 0). In the presence of 
the SOI at (/3 = 0, the last term is relevant for 5 > in 
addition to the first term. The first term is independent 
of T, whereas the last term increases with decreasing 
T as r2/f-2 j-Qj. ^ > 0. For Lp ^ 0, the second term 
is relevant for 5 < even in the absence of the SOI. 
Finally for 7^ and A 7^ 0, all terms contribute to 
the conductance. The argument above is summarized in 
Table H 

The total current through a Josephson junction is de- 
scribed by so called resistively and capacitively shunted 
junction model 



J 



C (fp 
Yeltt? 



1 



dif 



2eR{ip) dt 



+ Josin((^), (24) 



with C being the capacitance of a junction. In the present 
junction, the resistance R — \/u depends on ip. Thus 
Majorana fermion excitation would modify dynamics of 
a Josephson junction. The phase Lp may be determined 
self-consistently so that the current can be optimized. 
Such issue is a natural extension of this paper. 

In the case of opposite helicity configuration^ we also 
obtain the Josephson current and the conductance as fol- 
lows 



(J) = eA 



at 


2 




sin p 


TTV 



1 

K 



XlK 



(25) 



<7 sin^(y;/2) ^ „„,2 n 
= TT — hA^cos y-^J^e 



Gq 



K 



+ TrX3Kcos^{p^/2) + Xl sin^ [^)d^ 



T 
T 



2/K-2 



2K-2 



(26) 



For A = 0, the conductance proportional to sm'^{(p/2) is 
independent of temperature, which is in sharp contrast 
to that in the equal helicity case given in Eq. (|23|) . The 
behaviors of the conductance a are summarized in the 
Table D 



In summary, we have studied electric transport 
through a point contact which connects Majorana 
fermion modes appearing at the edges of two helical 
superconductors by taking into account the interaction 
among Majorana fermions and the spin-orbit interac- 
tion (SOI) at a point contact. By introducing a fic- 
titious fermion consisting of two Majorana fermions, 
the Majorana fermion model is transformed into the 
Tomonaga-Luttinger model. The application of the stan- 
dard bosonization technique enables to analyze low en- 
ergy physical phenomena of a Majorana fermion. It is 
found that several novel features appear due to the Ma- 
jorana fermions such as (i) the conductance is sensitive to 
the phase difference of two superconductors, (ii) tunnel- 
ing with SOI gives quite distinct behavior from that with- 
out SOI, (iii) the transport phenomena depend on rela- 
tive helicity of two superconductors as shown in Eqs. (j23p 
and (p6)l . and (iv) the interactions leads to the power- law 
temperature/ voltage dependences. 



[1] X. G. Wen, Phys. Rev. B 44, 5708 (1991). 

[2] C. Nayak, et. al, Rev. Mod. Phys. 80, 1083 (2008). 

[3] A. O. Gogolin, A. A. Nerseyan, and A. M. Tsvelik, 
Bosonization and Strongly Correlated Systems, Cam- 
bridge (1998). 

[4] D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001); 
S. Das Sarma, et.al, Phys. Rev. B 73, 220502(R) (2006). 

[5] G. E. Volovik, The Universe in a Helium Droplet, Claren- 
don Press, Oxford (2003); Y. Nagato, et. al, J. Phys. 
Soc. Jpn. 78, 123603 (2009); R. Roy: arXiv:0803.2868 
S. B. Chung and S.-C. Zhang, Phys. Rev. Lett. 103, 
235301 (2009). 

[6] S. Tewari, et. al, Phys. Rev. Lett. 98, 010506 (2007); T. 

Mizushima, et. al, Phys. Rev. Lett. 101, 150409 (2008). 
[7] L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 

(2008); L. Fu and C. L. Kane, Phys. Rev. Lett. 

102, 216403 (2009); A. R. Akhmerov, J. Nilsson, and 

C. W. J. Beenakker, Phys. Rev. Lett. 102, 216404 (2009). 
[8] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
[9] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 



(2005); B. A. Bernevig, and S. C. Zhang, Phys. Rev. Lett. 

96, 106802 (2006). 
[10] A. P. Schnyder, et. al, Phys. Rev. B 78, 195125 (2008). 
[11] X. L. Qi et. al, Phys. Rev. Lett. 102, 187001 (2009); 

M. Sato and S. Fujimoto, Phys. Rev. B 79, 094504 

(2009). 

[12] Y. Tanaka, et. al, Phys. Rev. B 79, 060505R (2009). 

[13] Y. S. Hor, et. al, Phys. Rev. Lett. 104, 057001 (2010), 
X. L. Qi, et. al, arXiv: 1003.5448 

[14] Y. Tanaka, et. al, Phys. Rev. Lett. 103, 107002 (2009); 
J. Linder, et. al, Phys. Rev. Lett. 104, 067001 (2010); 
K. T. Law, P. A. Lee, T. K. Ng, Phys. Rev. Lett. 103, 
237001 (2009); L. Santos, et. al arXiv:0 9l07592Tt J. Al- 
icea, Phys. Rev. B 81, 125318 (2010); J. D. Sau, et. 
al, Phys. Rev. Lett. 104, 040502 (2010); J. Nilsson and 
A. R. A khmerov, arXiv: 0912.4716; R. Shindou, et. al, 
larXiv: 1004.0750 

[15] J. C. Y. Tao and C. L. Kane, Phys. Rev. B 79, 235321 
(2009). 

[16] A. F. Ho and P. Coleman, Phys. Rev. Lett. 83, 1383 



5 



(1999). 

[17] C. Wu, et al, Phys. Rev. Lett. 96, 106401 (2006); 
C. Xu and J. E. Moore, Phys. Rev. B 73, 045322 (2006); 
Y. Tanaka and N. Nagaosa, Phys. Rev. Lett. 103, 166403 
(2009). 

[18] T. Giamarchi, Quantum Physics in one- dimension (Ox- 
ford Science PubUcations, 2004). 
[19] J. Solyom, Advances in Physcis, 28, 209 (1979). 



[20] H.-K. Kwon. et. al, Eur. Phys. J. B, 37, 349 (2004). 
[21] A. Furusaki and N. Nagaosa, Phys. Rev. B 47, 4631 
(2003). 

[22] C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 

1220 (1992); Phys. Rev. B 46, 15233 (1992). 
[23] M. P. A. Fisher and W. Zwerger, Phys. Rev. B 32, 6190 

(1985). 



